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Abstract. We are concerned with the convergence of a numerical scheme 
for the initial-boundary value problem associated to the Korteweg-de Vries- 
Kawahara equation (in short Kawahara equation), which is a transport equa- 
tion perturbed by dispersive terms of 3rd and 5th order. This equation appears 
in several fluid dynamics problems. It describes the evolution of small but finite 
amplitude long waves in various problems in fluid dynamics. We prove here 
the convergence of both semi-discrete as well as fully-discrete finite difference 
schemes for the Kawahara equation. Finally, the convergence is illustratred by 
several examples. 



1. Introduction 

This paper is concerned with the initial-boundary value problem of the Kawahara 
equation: 

(1.1) V>t ^^x ^xxx ~t~ ^xxxxx ? 

with initial condition 

(1.2) u(x,0) = f(x), for all x 
and the boundary condition 

(1.3) u(x,t) = u(x + l,t), for all x and t 

It is well known that the one-dimensional waves of small but finite amplitude in 
dispersive systems (e.g., the magneto-acoustic waves in plasmas, the shallow water 
waves, the lattice waves and so on) can be described by the Korteweg-de Vries 
(KdV in short) equation, given by 

(1.4) Ut = -uu x - u xxx , 

which admits either compressive or rarefactive steady solitary wave solution (by a 
solitary water wave, we mean a travelling wave solution of the water wave equations 
for which the free surface approaches a constant height as \x\ — > oo) according to 
the sign of the dispersion term (the third order derivative term). Under certain 
circumstances, however, it might happen that the coefficient of the third order 
derivative in the KdV equation becomes small or even zero. In that case one has 
to take account of the higher order effect of dispersion in order to balance the 
nonlinear effect. In such cases one may obtain a generalized nonlinear dispersive 
equation, known as Kawahara equation, which has a form of the KdV equation with 
an additional fifth order derivative term given by (1.1). The Kawahara equation is 
an important nonlinear dispersive equation. It describes solitary wave propagation 
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in media in which the first-order dispersion is anomalously small. A more specific 
physical background of this equation was introduced by Hunter and Scheurle [3], 
where they used it to describe the evolution of solitary waves in fluids in which 
the Bond number is less than but close to | and the Froude number is close to 1. 
In the literature this equation is also referred to as the fifth order KdV equation 
or singularly perturbed KdV equation. The fifth order term d^u is called the 
Kawahara term. There has been a great deal of work on solitary wave solutions 
of the Kawahara equation [4, 5, 8, 10, 11] over the past thirty years. It is found 
that, similarly to the KdV equation, the Kawahara equation also has solitary wave 
solutions which decay rapidly to zero as t — > oo, but unlike the KdV equation 
whose solitary wave solutions are non-oscillating, the solitary wave solutions of the 
Kawahara equation have oscillatory trails. This shows that the Kawahara equation 
is not only similar but also different from the KdV equation in the properties of 
solutions, like what happens between the formulations of this equation and the 
KdV equation. The strong physical background of the Kawahara equation and 
such similarities and differences between it and the KdV equation in both the form 
and the behavior of the solution render the mathematical treatment of this equation 
particularly interesting. The Cauchy problem given by (1.1), (1.2) has been studied 
by a few authors [6, 7, 9, 12, 13]. It has been shown that the problem (1.1), (1.2) 
has a local solution u G C([—T, T];H r (R)) if / G H r {R) and r > -1. This local 
result combined with the energy conservation law yields that (1.1) has a global 
solution u G C([— oo, oo]; L 2 (R)) if / G L 2 (R). Well-posedness results can be found 
in [6]. 

In this paper, we focus on the the derivation of convergent finite difference numer- 
ical methods for the initial-boundary value problem (1.1), (1-2), (1.3). The problem 
of analyzing convergent numerical schemes of course intimately connected with the 
mathematical properties of the Cauchy problem for the KdV equation, which is well 
developed in literature. As far as we are concerned, the first mathematical proof of 
the existence and uniqueness of the solutions of the KdV equation (1.4) was accom- 
plished by Sjoberg [2] in 1970, using a finite difference semi-discrete approximation. 
In [1], authors have considered a fully-discrete finite difference scheme for the KdV 
equation and showed the convergence of their scheme. In this paper, we are going 
to use a similar technique here. Note that Sjoberg's approach is based on a semi- 
discrete approximation where one discretizes the spatial variable, thereby reducing 
the equation to a system of ordinary differential equations. However, we also stress 
that for numerical computaions this set of ordinary differential equations should 
be further discretized in order to be solved. Thus in order to have a completely 
satisfactory numerical method, one seeks a fully discrete scheme that reduces the 
actual computation to a solution of a finite set of algcbric equations. In the present 
paper, we consider both semi-discrete and fully-discrete cases to prove the following 
theorem: 

Theorem 1.1. If f(x) is a 1 -periodic function and fifth derivative of fix) belongs 
to L 2 (R), then there exists a unique solution of the problem (1.1), (1.2) and (1.3), 
i.e., there exists a function u G H 5 (R) with Ut G L 2 {R) for which the euality (1.1) 
holds as an L 2 equality. 

The numerical computaion of solutions of the Kawahara equation is rather capri- 
cious. Two competing equations are involved, namely nonlinear convective term 
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uu x , which in the context of the equation ut — uu x yields an infinite gradient in fi- 
nite time even for smooth data, and the linear dispersive terms u xxx , u xxxxx , which 
in the context of the equation ut = u xxx + u xxxxx produces hard to compute dis- 
persive waves, and these two effects combined makes it difficult to obtain accurate 
and fast numerical methods. Most of the finite difference schemes will consist of 
a sum of two terms, one discretizing the convective term and one discretizing the 
dispersive terms. These two effects will have to balance each other, as it is known 
that the Kawahara equation itself keeps the Sobolev norm H s (s > —1) bounded. 

The rest of the paper is organized as follows: In section 2, we consider the semi- 
discrete scheme for the initial-boundary value problem corresponding to (1 .1). At 
first, we state some of the well-known Sobolev type estimates and then we showed 
the convergence of the semi-discrete scheme. Both local and global existence has 
been proved. In section 3, we consider a fully-discrete semi-implicit scheme for the 
initial-boundary value problem corresponding to (1.1). We have used explicit dis- 
cretization for the "nonlinear" term and implicit discretizations for the "dispersive" 
terms. Convergence of the fully-discrete scheme has been shown in this section. In 
section 4, we have shown the uniqueness of the solution to the initial-boundary 
value problem given by (1.1). Finally, in section 5, we have justified the conver- 
gence results by several numerical experiments, namely Soliton experiments. We 
also compare our scheme, based on the fully-discrete semi-implicit scheme, with the 
existing results in literature. 

2. Semi Discrete Analysis 

2.0.1. Local existence: Let Ax be a small positive number and define a grid on 
the X-axis to be the set of gridpoints Xi = iAx for i £ Z. We will denote the 
value of our difference approximation at — Ui. Since we are in the 

periodic case, for simplicity we assume a unit period and that 1/Ax £ N. In that 
case Ui+N — ^ for j g Z. To simplify the notations, we will introduce the finite 
difference operators: 

D-Ui = *-r— (Ui - D + Ui = -r— (U i+1 - Ui), D Ui = — r— - Ui-l). 

Ax Ax 2Ax 

We will also use the following notations: 



Il/H 2 = (/,/) and (/,<?)= / f(x)g(x)dx, 



(i 



and in the space of gridfunctions ( a discrete, possibly complex valued, function 
defined on the grid ), we define the scalar product and the norm by 



(f,g)h = hJ2f(x l )g(x l ) and ||/|j£ = (/,/)/,. 



Later we need the following Lemmas, proven in [2]. We begin with a well-known 
Sobolev-type Lemma: 

Lemma 2.1. Let a and r be integers such that < t < a. Then for every 
constant e > there exists a constant c(e) such that for all functions y, sufficiently 
differ entiable on < x < 1, 



(2-1) max& 2 <e 
o<x<i , dx T ~ 



dx° 



2 



c(e)\\y\\ 2 
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and 
(2.2) 



d T y 


i 

< e 


d°y 


dx T 




dx a 



c(e)\\y\\ 



Lemma 2.2. Let T\ and T2 be nonnegative integers with t\ + T2 = r and ip be a 
function of the form 



Then 





d T Tp 




dx T 



< \\D^D?^\\ = \\D^D^ilj\\ h < 



dx T 



Lemma 2.3. Let a and r be integers such that < r < a , and z be any gridfunc- 
tion. Then for every constant e > there exists a constant c(e) independent of h 
and z such that, 



(2.3) 



\D T + z\\ 2 h <e\\D° + z\\ 2 h +c{e)\\z\\l. 



Remark 2.1. The inequality (2.3) can be modified as follows: One can replace 
by any operator D T + D^ 2 D I 3 with r = t\ + T9 + T3 . The right-hand side can 
7? 



likewise be replaced by D^ 1 D°i 2 with a = a\ + 02 . However, it can be proven by 



counterexamples that the lemma is not true if the right member contains Do ■ 
Remark 2.2. From the above Lemmas the folloing result follows immediately 
(2.4) max|A)^| 2 < e\\D%D 2 _zf + c{e)\\z\\l, 



where e and c(e) have the same properties as before. 

Now we are in a position to state the semi-discrete scheme of the Kawahara 
equation (1.1), given by 

(2.5) {u l ) t ^~ l -[u l D u l + D uj]-D^D\u l + DlD 2 _u l , i = 1, 2, N, 
with the initial condition 

(2.6) Ui(0)=f( Xi ), » = 1,2,...,JV 
and the boundary condition 

(2.7) Ui(t) = Ui + m(t) for all i and t. 

We will first prove the local existence of a solution for t > 0. The case t < will 
be treated later. 

Theorem 2.1. There exists a time T\ > and constants fc,, i — 0, 1,2,3,4 inde- 
pendent of Ax but dependent on f{x) and its derivatives of order five and lower, 
such that 



(2.8) 
(2.9) 
(2.10) 
(2.11) 



IK-,*)llh < fco, f° r al1 1 

\u(x u t)\ < h, 0<t<T u 

\D_D\u{-,t)\\ h <k2, 0<t<Ti 

\D%D 2 _u(;t)\\ h <kz, 0<t<Tx 



for all i 



KAWAHARA EQUATION 



5 



and finally with, 

du{x,t) 

(2.12) K-.*)llh<*4. 0<t<Ti. 

Proof. Multiply the above equation (2.5) by AxUi and summing over all i, we have 

{u,u t )h = --[(u,uD u) h + (u,D u 2 ) h } - (u,D-D\u) h + (u,D\D 2 _u) h 

Aj; „„ „ „2 Ax I, „ „o 1 1 2 

= D+D-U t \\D_D 2 u 111 

2 h 2 + lift 

where we are using the following identities: 

{u,D_D 2 + u) h = ^{u,D_D 2 + u) h - l -{u,D+D 2 _u) h 

1 
2 

= ^\\D + D_uf h , 

since {u, D+D^iijh — —(UjD-D+ujh, because {u,D + u)h — —(u,D-u)h, in the 
first line, and 

£>+-£>_ = AxD_D+ = AxD+D_. 
In a similar manner, we find that 

{u,D\D 2 _u) h = l -{u,DlD 2 _u) h - ~(u,D 3 _D 2 + u) h 
= ^(u,D 2 D 2 + (D + -D.)u) h 

Ax I, 2 || 2 

= — \\D-D±U , . 
O II + lift 



(2.13) = -(«, D+D_ (D+-D-)u) h 



So, finally we conclude that 



\H;t)\\l + ~J o (\\D + D_ U (r)\\l+\\D-D 2 + u(T)\Q dr 

<h(-,0)||^||/|| 2 , <2 / f 2 (x)dx = k 2 , 
Jo 

from which (2.8) follows. 

Keeping in mind that v — ut, from the equation (2.5) using the triangle inequality 
and Lemma 2.3, we get the following inequality 

\\D\D 2 _u\\ h < \\v\\ h + i || [uD u + D Q u 2 ]\\ h + \\D_D\u\\ h 

< Uplift + max\D u\ \\u\\ h + (e \\D 3 + Dlu\\ h + c(e) \\u\\ h ) 



< 



\\v\\ h +\\ U \\ h (e\\D d + Diu\\ h + c(e)\\u\\ h ) 



+ (e\\D d + Diu\\ h +c(e)\\u\\ h ) 
Now we can use the fact that \\u\\ h < fco, and choose e such that 
(2-14) \\D 3 + D 2 _u\\ h < Vl \\v\\ h + V2 , 
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where v\ and v 2 are constants independent of Ax. To get a bound for \\v\\ h , we 
will proceed as follows. After differentiating both sides of the equation (2.5) with 
respect to t, we get the following equation 

(2.15) (vi) t = -g[«jD «i + UiDoVi + 2D {u lVl )\ - D^D 2 +Vl + D\D 2 _ Vl . 

Now multiplying equation (3.15) by hvi and summing over i, we have 
1, 



(v,v t )h 



< 



[(v, vD u) h + (v, uD v) h + 2(v, D uv) h ] 
- {v,D-D 2 + v) h + {v,D%D 2 _v) h 
(v 2 ,D u) h + (v, D uv) h \ - (v, Drfv) h + (v, D%D 2 _v) h 



1 



A 



1 



< -\(v ,D u)h + — (v,D + uD v) h + -(v,Vj-iD + u) h \ 



3 

< C 



max \D u\ \\v\\ h + max \D+u\ \\v\\ h 



j t \ni<v 3 \\vf h +u4,\\v\\i, 



<C\\v\\ 2 h [e\\DlD 2 _u\\ h + c(e)k Q ], 
where C is a constant independent of Ax. Hence we can conclude that, using (2.14) 

(2.16) 

where and are constants independent of Ax. So, from (2.16), it is clear that 

0<\\v\\ 2 h <y, 
where y is the solution of the initial-value problem 

(2-17) ^ = w*+w, y(o) = m>\\v(-M\l- 

It can be shown that the solution of (2.17) is finite for 

t < t 



-log(l + ^). 



There is no possibility of choosing constants scuh that — oo. However, we 
can always find a constant C such that (2.12) holds for < t < T% = too/2. The 
estimate (2.11) then follows from (2.14). Also the estimates (2.9) and (2.10) follows 
from Lemma 2.3 and (2.11). This completes the proof of the theorem. □ 

Now we can use the Theorem 2.1 to conclude that the Kawahara equation (1.1) 
has a solution u(x,t) in < t < T\. To do that, let us first assume that u^ x is 
a grid function, i.e., a function uax '■ [0,T] — > R N , where u& x {xj,,t) = Ui(t), such 
that m = un- We will also denote ipAx(x,t) to be the Fourier series of UAx(x,t). 
Now we are ready to use the Lemma 2.2 and 2.1 to conclude immediately that 

2 

< Cl, 



(2.18) 
(2.19) 
(2.20) 



&>pAz(;t) 



dt 

dx 3 

5 1pAx(;t) 



dx 5 



<c 2 , 



<c 3 , 



0<t<T u 
< t < Ti, 
Q < t < T 1} 
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where ci, C2 and C3 are constants depending on f{x) but not on Ax. Now the above 
inequalities (2.18) and (2.19) imply that the sequence Wa^Ia^x) ls bounded and 
equicontinuous in both x and t. The Arzela-Ascoli theorem guarantees the existence 
of a subsequence of {ip&x}& x> o which converges to some function u(x, i) as Ax —> 0. 

The rest of the argument is standard. Both 9 ^Jt * and d q x $ x converges to and 

respectively in the 1? sense. It is also easy to see that u{x,t) is a solution of 
the Kawahara equation because of the very definition of ipAx(x,t). 

2.0.2. Proof of Global existence. In this section we are going to prove the existence 
of the global solution of the equation (1.1). At first we are going to prove the 
following Lemmas, which we will use later. 

Lemma 2.4. Let u(x,t) be a solution of the problem (1.1). Then there exist a 
constants ai,ct2 such that 

(2.21) ( u 2 {x,t)dx= f u 2 (x,Q)dx= f f 2 dx = ax 
Jo Jo Jo 

(2.22) jf * - u\ - .L) dx = £ Q/ 3 - f - Z" 2 ) dx = a 2 

Proof. To prove (2.21) we start by multyplying the equation (1.1) by u and integrate 
by parts in space, yields 

tJ/Uj-f (jLOb / IXi < \Jb^ ' UL'JJj 'jq I ^-^^-^ oo oo oo oo oo 

1 1 f 1 1 f 1 

(-u 3 ) x dx- / (uu xx - -u 2 x ) x dx - I u x u xxxx dx 
& Jo z Jo 

~J (^u 3 ) x dx-J (uu xx - ^u 2 x ) x dx 

1 1 2 

{u x 1l xxx ~^^xx)x dx 0. 

Since all the boundary terms vanishes due to the periodic boundary condition. 
Hence this established (2.21). 

To prove (2.22), we start by multyplying (1.1) by u 2 and integrate by parts in 
space, yields 

1 2 f 1 3 2 2 

11 Jo 

1 1 f 1 f 1 

(-u i ) x dx + 2 (uu x )u xx dx-2 {uu x )u xxxx dx 

4 Jo Jo 

2 / [ 1H V'xxx ~t~ ^xxxxxl^xx dx 

Jo 

2 / [ Ut n xxx -\- u X xxxx\^xxxx dx 
Jo 

,1 ,1 ,1 

2 / r ^ J t,02 ^^33 ~ 2 / ^'OOOO O0XO0 I ^ / ^^OOOO ^'OOOOOOOO 00 

Jo Jo Jo 
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f 1 f 1 

2 / Utx^xxx dx -\- 2 / UxxxUxxxx 

JO JO 

= 2 / u tx u x dx + 2 u txx Uxx dx. 
Jo JO 

From this we can conclude that 

- 

dtJo L 3 

and consequently (2.22) follows from the above equation 



dx 2 / Uxxxx^xxxxx 




dx 



a - u x — u xx ] dx = 0, 



□ 



Lemma 2.5. Let u(x,t) be a solution of the problem (1.1). Then there exists a 
constant a such that 



(2.23) 
(2.24) 



max \ u x (x, t)\ < a 



< e 7t 



-// -/ +/ 



du 

at' 



Proof. From (2.21), it follows that 



< 



3 Jo 



\u 3 \dx+\\u x \\ +\a 2 \ 



< 3 ( c ( e ) IN + e \\ u xx\\) \\u\\ 2 + (c(e) \\u\\ 2 + e \\u xx \\ 2 ^ + \a 2 \ 

= 3 ( c ( e V«i + e \\ u xx\\) Oil + (c(e)ai + e ||w x:c || 2 ) + |a 2 |- 
Now we can rewrite the above inequality in the following form 
(2.25) a\\u X x\\ 2 ~ b\\u X x\\ ~ c <0, 

for some constants a, b, c, where a — 1— e, b = \cn\e and c = |c(e)y / a7+c(e)o;i + |a2|. 
Now it is easy to see that (2.25) gives, 

2 ,2 



'a - 



2-y/a 



< c 



4a' 



From the above relation, it is clear that < 0:3, for some constant a.3. Again 

using the interpolation inequality, \\u x \\ < (c(e) ||u|| + e Hunll), we can conclude 
that \\u x \\ < «4, for some constant 04. Also we can use a similar type interpolation 
inequality to conclude that (2.23) holds. 
Now the function v(x, t) 



tj| satisfies, 



dv 



Multiplying the above equation by v and integrating in space yields, 

\a\ " W " 2 = VUx ^ ~ ( V,UVx ^ ~ ( V ' Vxxx ) + ( v > v xxxxx) 



Now (2.23) gives, 



|||t,|| 2 <maxK|||t,|| 2 <C||t-|| 2 , 
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which implies 



n 2 
v\ = 



du(-,t) 



dt 



du(-,0) 



dt 



alt 



-ff' - f + f 



i.e. we established (2.24) 
Now to get a bound on \\u 



a 



1 d 

2~dt 



u xxx ^ 



, , we proceed as follows: 
l 

dx 



(i 



' 'XXXX a xxl 



dx = 



u x t dx 



Uxt{u t 



vv T dx 



(-v 2 ) x dx + 



~) dx 



uu x u x t dx 



uu x u xt dx 



= / vu x dx + vuu xx dx + /3 



1 d 
2dt 

1 d 
~ 2dt 
1 d 
'2~dt 



. dx 



< (maxl^Ddlu 

1 d 



+ 



2dt 



\v\\) + (maxluDdlu^H ||u|| 



. dx 



Since the use of the interpolation inequality allows us to conclude that all the terms 
max max \ u x \, \\u x \\ , Hit^x II are bounded by some constant, hence from the above 
relation we end up with 



d_ 

dt jo 

Therefore we conclude that, 
(2.26) \\u xxx \\ 



1 



1 



u 2 xx dx < a\\v\\ 



< Cie 7t ( 



-//' - / + /" 



Again we can now use the Kawahara equation (1.1) and triangular inequality to 
conclude that From the above relation we conclude that, 



(2.27) 



|| Uz 



< de 7 *( 



-//' - / + /" 



C 2 



where C\ and C2 are constants. One can see that the bound (2.24) guarantees that 
|jjr is square integrable for every t and (2.26), (2.27) that the problem (1.1) with 
the initial function u(x,T 1 ) instead of f(x) has a solution for T\ < t < T 2 = 2Ti. 
Consequently, we get a solution of the Kawahara equation for < t < T. Now to 
obtain the existence of sloutions for all t > 0, we will repeat the extension procedure. 
For this purpose suppose that existence can be proven only in < T < 00. Now if 
we look at the expression for t^, we find that only yo depends on t. But yo can, 
because of (2.24), be chosen to hold in the whole interval < t < T. Consequently, 
if we consider problem (1.1) with f(x) = u(x,t) for some r sufficiently close to 
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T, we can, by using the local procedure, get existence for values of t lying outside 
< t < T. 

Remark 2.3. To prove the existence in the lower half-plane, we proceed as follows: 
Consider the following equation 

(2.28) — UU X -\- U X xx ^XXXXX 1 

which arises when making the transformations x — > (1 — x) and t — > — t. Now if 
we substitute (2.28) for (1.1), then we can prove existence of solution in the lower 
half plane. This is possible since we get an interval Ti < t ; < in Theorem 2.1. 
Consequently, this concludes the proof of the existence part of Theorem 1.1. 

3. Fully Discrete semi-Implicit Scheme: 

We know that in order to have a completely satisfactory numerical method, one 
must seek a fully discrete scheme that reduces the actual computation to a solution 
of a finite set of algcbric equations. So, here we will consider a semi-implicit fully 
discrete scheme and show the convergence of the solution of that fully discrete 
scheme to the solution of (1.1). To do that, let us select a time step At > 0, and 
write t n = nAt. The value of our differnce approximation at (j Ax, nAt) will be 
denoted by w™ in the fully discrete case. To simplify the notation, we introduce the 
finite difference operator: 

We propose the following semi- implicit fully-discrete aproximation to (1.1), given 

by 



a" : ~ 1 = a 1 ; - A/ v; /;.,«; + A>K) 2 ] - Aii) /)-//;; •' + aid: d 2 



3 3 Q 

(3.1) 3 | 
where Uj = -(uj+i +itj_i). 

Now keeping in mind that DqUj = lu^D^Uj, we can rewrite the above scheme as 

(3.2) u] +1 = u'j - Atu]D u] - AtD_D 2 + u n + 1 + AtD% D 2 _ u] +1 . 
with the initial condition 

(3.3) u? = /(si), i = l,2,...,JV, 
and the boundary condition 

(3.4) uf = u™ +N for all i and n. 

3.1. Local existence of solution: To show the local existence of solution to (1.1), 
we will use similar arguments to the ones used in the semi discrete case. First we 
will state the main theorem, similar to Theorem 2.1. 

Theorem 3.1. There exists a time T > and constants hi, i — 0,1,2,3,4 inde- 
pendent of Ax but dependent on f(x) and its derivatives of order five and lower, 
such that 

(3.5) \\u n \\ h <k , 0<nAt<T 

(3.6) \u n (xi)\<ki, 0<nAt<T 
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(3.7) \\D_D 2 + u n \\ h <k2, 0<nAt<T 

(3.8) \\D 3 + D 2 _u n \\ h <k 3 , 0<nAt<T 
and finally with, 

v n = D+u n -\ neN 

(3.9) \\v n \\ h <h, < nAt < T. 

Before giving a proof of the above theorem, we will first prove the following usefull 
Lemma: 

Lemma 3.1. Let u™ be a solution of the difference scheme (3.2) and define v n = 

D^u n ~ x . Then the following two estimates hold 

(3.10) 

«+i|| 2 _i_ a + a-». Ml r> r> 2 „«+i|| 2 . II n . n „™+ 1 ll 2 _i_ _ II r>„„,™ll 2 I <? ll,, n ll 2 



p || + At Ax (J|L>_L» 2 u n+i || + ||D + £>_u n+i || + — \\D u n \\" ) < \\u n \\ 
also we have, 

\\v n+1 \\ z + AtAx [ \\D+D-v n+1 \\' + \\D-Dl v n+1 \\" + \\Dv 



< \\v n \\ 2 +CAtm&x\Du n \ \\v n \\ 2 . 



(3.11) " M V M ' " " " S A 

provided the CFL condition 

<W ^Ml(i + ^K||)<|, 

Proof. First of all, just to avoid clumsy notations, we will drop the index j from 
our notation, and use the notation it, Du for Uj, DgUj respectively where j is fixed. 
We first study the "Burgers" term AtuDu. Let u be a gridfunction and set 

(3.13) w = u-AtuDu. 

Set A = At/ Ax, and we will use the following CFL condition: 

/o 1 a\ At ., At „ At „ . At ... 3 

(3.14) _ max |u|(l + 24— max |„|) < \\u\\ (1 + 24^ \\u\\) < - 

Multiplying (3.13) by u, we have 

-w 2 = l u 2 - At^uDu 2 + -(w- u) 2 
2 2 2 2 y ' 

= l u 2 - At-uDu 2 + At 2 u 2 (Du) 2 + -(u 2 - u 2 ). 
2 2 2 

We will use the following relations 

1 l{a + b)(a 2 -b 2 )= 1 -(a 3 -b*)-±(a-b)\ 



and 



\ {a + b) 2 -\{a 2 + b 2 ) = -\{a~b) 2 . 



For a gridfunction, these implies 



uDu 2 = 1 Du 3 - ^^(Du) 3 , 
3 12 V ' 

u 2 -u 2 = -—(Duf. 
4 
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Therefore, 

1 , 1 -„ At o . 9 n9 AtAx 2 , x o Ax 2 , x9 

1 - At At 2 

(3.15) = * u 2-^ Du S-^L {Du f 

+ —(Duf ( XAxDu + 2A\ 2 u 2 - - 
24 V 2 / 

Now summing (3.15) over j and using the CFL condition (3.14) we conclude, 

(3.16) l^f + ^n^f < M a. 

o 

Now we are in a postion to study the full difference scheme by adding the "Airy 
term" AtD-D\u^ +l , and the "Kawahara term" AtD\D 2 _v!] +1 . Thus the full 
difference scheme (3.2) can be written as 

(3.17) v = w - AtD_D\ v + AtD\D_ 2 v. 
Multiplying the above equation (3.17) by v, and summing over j, we have 

\\M 2 + J2 - w ^) 2 = § ll^ll 2 — At(«, D.D\v) + At(v, D\D 2 _v) 

j 

1 „ ,,2 AxAt „ „ „ „ 2 AxAi 
= 2 INI --j 

Hence finally we have, 



HI — H-D+^-f H Fi — 11-°-°+^ 



(3.18) ||u|| 2 + AxAt \\D+D-v\\ 2 + AxAt \\D-D 2 + v\\ < \\w\\ 2 

Now combining two results coming from (3.16) and (3.18), we conclude that 
(3.19) 

|U" +1 f + AxAt\\D + D-u n+1 \\ 2 + AxAt\\D + D 2 _u n+1 \\ 2 + ^\\Du n \\ 2 < \\u n \\ 2 

8 

which is exactly (3.10). To prove (3.11), we will proceed as follows: 
First we define 

v n = D+u n -\ ueffo 

Now the gridfunction satisfies 
(3.20) 

v n+l = yn_ At ^n Du n + gn^") + A ^w"D W " - AtD_D 2 + V n+1 + AtD 3 + D 2 _ . 

Set 

At 2 

w = v- AtD(uv) + —^Dv 2 

We proceed as before and multiply this by v to find 

1 2 1 -„ At 2 / . , At , x 2 
-w 2 = - v * + — - ~2 

Ax 2 

- At{uvDv + v 2 Du) + At 2 v 2 Dv ^~{Dv) 2 

8 
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Now we will use the following relations 



D(uv) - ^Dv 2 ) < (D(uv)) 2 + At 2 v 2 (Dvf 



< 2u 2 (Dv) 2 + 2v 2 (Du) 2 + M 2 v\Dv), 

v 2 Dv = -Dv 3 - ^—(Dv) 3 , 
3 12 y ' ' 

1 1 Ax 2 

uvDv + v 2 Du = -D(uv 2 ) + -v 2 Du (Dv) 2 Du. 

2 2 8 

Using this 

-u; 2 = -v 2 - AiD -ot 2 v 3 v 2 Du (Dv) Du 

2 2 \2 3 J 2 8 1 ' 

(3.21) + At 2 [2u 2 {Dvf + 2v 2 (Du) 2 + At 2 v 2 (Dv) 2 - ^{Dvf^j 

At 2 

Now our aim should be to balance all the positive terms with Ax 2 (Dv) 2 . For that 
we will proceed as follows: 

Ax 2 (Dv) 2 < 4v 2 , 
At 2 v 2 (Duf < CAtv 2 \Du\, 

At 2 v 2 < 2 ^max|u™| 2 

Ax 2 \Dv\ < -\u\, 
A 

where the constant C depends on ||u||. Now using the above inequalities in (3.21), 
we get 

Z „.2 A + n I „.„.^ „.3 



maxlu"- 1 ! 2 



-w 



= -v 2 - AtD -uv 2 - —v 3 + CAtv 2 \Du\ 



2 2 

+ X 2 Ax 2 (Dv) 2 (4max|u™| 2 + 2max|u n - 1 | 2 + -^-max\u r ' 
8 v ' 

Now we choose At in such a way that, 

(3.22) A 2 (4max|u"| 2 + 2max|u n - 1 | 2 ) + ~max|u™| < ~. 
Now after summing over j, we have 

At 2 

(3.23) || ^ || 2 + — ||^|| 2 < ||v|| 2 + CAtm^\Du\ \\vf . 



Now using that 



v n+1 =w- AtD+D 2 _v n+1 + AtD 3 ,D 2 _v n+1 , 
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we get 
(3.24) 



\v n+1 f + AxAt \\D+D-v n+1 f + AxAt \\D^D\ v n+1 f + — \\Dv n \\ 2 

8 

< ||w"|| 2 + CAtmax|L>u Tl | ||u™|| 2 . 



Which proves the second part (3.11) of the Lemma. □ 
Now we are ready to prove Theorem 3.1. 
Proof, (of Theorem 3.1) 

First note that, from the definition of v n , (3.2) can be rewritten as, 

(3.25) v n+1 = D+u n = ^D + D_u n - u n Du n - D_D\u n+1 + D 3 + D 2 _u n+1 

2 A 

Therefore, 

At 

\\D 3 + D 2 _u n+1 \\ < \\v n+1 \\ + ||u n Du"|| + -f \\D + D-U n \\ + \\D_D 2 + u n+1 II 

< ||u n+1 || +max\Du n \ \\u n \\ + \\D_D 2 + u n+1 \\ + C 

< \\v n+1 \\ +C*i \\D\D 2 _u n+1 \\ + C 2 \\D-D 2 + u n \\ 

< \\v n+1 \\ +Ct \\D 3 + D 2 _u n+1 \\ 

+ (AtC 2 \\D_D 2 + v n+1 \\ +C 3 \\D-D 2 + u n+1 \\) 

< \\v n+1 \\ +Ci \\DlD 2 _u n+1 \\ 

1 

+ (jK|| 2 (l + C*Atmax|L>u"|))\ 

where we have used (3.24) to estimate At ||Z?_.D+?; ra+1 1| . Now by the CFL condition 
(3.12), Atm&x\Du\ < C for some constant C. Hence, finally we have 

(3.26) \\D 3 + D 2 _u n+1 \\ <c + Cl \\v n+1 \\+c 2 \\v n \\, 

for some constants cq, c±, c 2 that are independent of Ax 
Now using this result in (3.24), we have 

IK +1 f < ikii 2 + At (d x \\v n \\ 2 + d 2 (\\v n \\ 3 + \\v n f ll^- 1 !))) , 

for constants d\ and d 2 . Set a n = ||w"|| 2 , so that 

a n +i < a n + At (d\a n + d 2 (a 7 l + a n o^_i 



Now let a be the solution of the differential equation 

da 3 

— = d\a + d 2 a 2 , «o = a > 0. 
dt 

The solution has a blow up time 

roo =<° + ^( i+ | 

Furthermore, for t < T° 
nAt < T°°(a, 0), a n < a(nAt; (0, a)). This clearly holds for n = and n = 1 (since 
a is increasing in t). Assuming that the claim holds for integers up to n, we get 

a n +i < a(nAt; (0,a)) 
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+ Atd x a{nAt; (0,a)) 

+ Atd 2 (a(nAt; (0,a))s + a(nAt; (0,a))a((n - i)At; (0,a))^ 

< a{nAt; (0,a)) + At (dta(nAt; (0,a)) + 2d 2 a(nAt; (0,a))^ 
<a((n + l)At; (0,a)). 

Hence for t <T = T°°/2, < C for some constant C independent of Ax. From 
this all the results in Theorem 3.1 follows. □ 

Therefore, we can follow exactly the approach of the semi discrete case to con- 
clude that uax converges in L? to some function u(x,t) for t < T. Furthermore, 
we have that 

\\ut\\ < ci, Hua^zll < c 2 and , 

Uxxxxx 1 1 — C3 

Therefore we can show that u(x,t) satisfies (1.1). 



v xxxxx 



4. UNIQUENESS 

To prove the uniqueness, let us assume that there exists two solutions u(x,t) 
and v(x, t) of the problem (1.1). Then the function w — u~v satisfies the following 
equation 

w t = -{uu x - vv x ) ~ w xxx + w x 
(4.1) = -wu x - vw x - w xxx + w xxxx 

w(x,0) = 0. 

Hence, by taking inner product of the above equation with w, we have 
. . 19||w|| 2 . . . . , . . . 

(W,W t ) = ^ = -(W,WU X ) - (W,VW X ) - {W,W XXX ) + {W,Wxxxxx) 

= -(w 2 ,u x ) + (w 2 , v x )/2, 
by periodicity. Now use of estimate (2.23) implies, 

<CNI 2 , 

for some constant C. Now as ||u;(-,0)|| 2 = 0, it is clear that ||iu(-,t)|| = for all t 
which consequently implies uniqueness. 



5. Numerical experiments 

5.1. Numerical experiment 1. The fully-discrete scheme given by (3.2) have 
been tested on a suitable numerical experiment in order to demonstrate its effec- 
tiveness. It is well known that a soliton is a self-reinforcing solitary wave that 
maintains its shape while it travels at a constant speed. Solitons are caused by a 
cancellation of nonlinear and dispersive effects in the medium. Several authors [ 
see [14], [15], [1] ] have studied the soliton experiments in the context of both KdV 
and Kawahara equation. Here we are interested in the soliton experiment for the 
Kawahara equation only. We are going to compare our result with the results given 
by [14]. 
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Although all the numerical experiments performed in [14] based on the equation 
given by 

(5.1) V>t ~t~ ^x ~t~ 1LU X -\- U X xx — ^xxxxx i 

but it is indeed very easy to see that (1.1) and (5.1) are completely equivalent by 
way of simple change of variables. In [14], authors have considered the following 
scheme for the Kawahara equation (1.1) 



- u? 1 



D\u?] 2 + Au? +1 = 



At 2 

where A = D + D + D_ — D + D + D + D_D_. In the case of Kawahara equation given 
by (1-1), if we consider the initial function in the following form given by 

then it is known from [16] that the explicit solution is given by the following trav- 
elling wave 

105 , 4 / i , m 



u(x.t) = — — sech — 7^(x — —— — c) 
169 V2yi3 169 ; 

This result can be verified through substitution. 

Since we know that the behaviour of the exact solution for Kawahara equation, 

mainly which remains its shape as time grows, it will be interesting to see how the 

numerical solution given by the scheme (3.2) evolves with time. We will use the 

following notations: UK scheme - scheme described in this paper and JMO scheme 

i 

- scheme described as in [14] and |M|; 2 = (Aa; J^k U D 2 ■ ^ n order to compare 
with the existing scheme given by [14] , we present the I 2 errors on a computational 
domain [—40,40], between exact solution and the solution generated by the UK 
and JMO schemes in table 5.1. 



Mesh points 


UK 


JMO 


4000 


2.7e-3 


1.2e-3 


8000 


1.4e-3 


6.2e-4 


12000 


9.2e-4 


4.2e-4 


16000 


7.0e-4 


3.0e-4 



Table 5.1. Numerical Experiment 1: I errors between exact and 
simulated solutions at time t = 10 for both U K and JMO schemes 



In the following figures we show the behaviour of the numerical solutions at 
different times. In this case we have used a domain [—20, 50], 5000 mesh points and 
a CFL number 0.75. We will compare our results with the results given by [14]. 

5.2. Numerical experiment 2. Now we will move on to the second numerical 
experiment. Here instead of one soliton we have considered two solitons, given by 
the initial data 

«(*, 0) = lg {sec* (^(* - 20)) + \sech* (-^(* - 60)) } . 
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: Exact, : UK scheme : JMO scheme. : Exact, ; UK scheme : JMO scheme 



Figure 5.1. Top Left: Exact and Numerical solution at time 
t = 30 for both UK scheme and JMO scheme; Top Right: Exact 
and Numerical solution at time t = 60 for both UK scheme and 
JMO scheme; Bottom Left: Exact and Numerical solution at time 
t = 90 for both UK scheme and JMO scheme; Bottom Right: Exact 
and Numerical solution at time t = 120 for both UK scheme and 
JMO scheme. 



This will essentially corresponds to the superposition of two solitons with differ- 
ent speeds given by the nonlinear term uu x of equation (1.1). We have used 10, 000 
points in space in the interval [—100,100]. Also we have run both schemes upto 
time t — 50. The following plot shows the behavior of the numerical solutions, 
where we see oscillatory structures of solitons. 




Figure 5.2. Left: Numerical solution for UK scheme; Right: 
Numerical solution for JMO scheme. 
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